function [ data ] = MockWHGaussWidinessByEnergy( lam, E )
%MOCKWHGAUSSWIDINESSBYENERGY Summary of this function goes here
%   Detailed explanation goes here
    
    Aampl = 1.5;
    Acenter = 530;
    Asigma = 70;
    Ae = sqrt(E);
    Agauss = GaussE(lam, Aampl, Acenter, Asigma, E)' ;
    Aspektras = Agauss;


    Bampl = 1.5;
    Bcenter = 400;
    Bsigma = 40;
    Bspektras = GaussE(lam, Bampl, Bcenter, Bsigma, E)' ;


    Campl = 1.5;
    Ccenter = 300;
    Csigma = 60;
    Cspektras = GaussE(lam, Campl, Ccenter, Csigma, E)' ;

    data = Aspektras + Bspektras + Cspektras;

end

function [ output ] = GaussE( x, ampl, center, sigma, E )
 sigmaBot =  2 * (E * 10).^ 2;
 A = -((x - center).^2); 
 eX = size(E, 2);
 output = repmat(A, eX, 1);
 sigmaBot = repmat(sigmaBot, size(x, 2), 1)';
 output = output ./ sigmaBot;
 
 output = ampl*exp(output);
end


